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Ground state energies and wave functions of quartic and pure quartic oscillators are calculated by 
first casting the Schrodinger equation into a nonlinear Riccati form and then solving that nonlinear 
equation analytically in the first iteration of the quasilinearization method (QLM). In the QLM 
the nonlinear differential equation is solved by approximating the nonlinear terms by a sequence of 
linear expressions. The QLM is iterative but not perturbative and gives stable solutions to nonlinear 
problems without depending on the existence of a smallness parameter. Our explicit analytic results 
are then compared with exact numerical and also with WKB solutions and it is found that our ground 
state wave functions, using a range of small to large coupling constants, yield a precision of between 
0.1 and 1 percent and are more accurate than WKB solutions by two to three orders of magnitude. 
In addition, our QLM wave functions are devoid of unphysical turning point singularities and thus 
allow one to make analytical estimates of how variation of the oscillator parameters affects physical 
systems that can be described by the quartic and pure quartic oscillators. 

PACS numbers: 03.65.Ca, 03.65.Ge, 03.65.Sq 



I. INTRODUCTION 

A basic nonrelativistic quantum mechanics problem is 
to solve the Schrodinger equation with a potential V{x) 
that governs motion of a given physical system. The first 
two terms of the power expansion of a one-dimensional, 
even potential around an equilibrium position arc 
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where x is the deviation from an equilibrium position. 
The above potential describes the dynamics of a great 
many systems that deviate from the idealized picture of 
pure harmonic motion. When both g and A are nonzero, 
we call this potential a "quartic" or quartic anharmonic 
oscillator; whereas, when g — Q with nonzero A it is 
dubbed a "pure quartic" oscillator. In addition to pro- 
viding an excellent description of spectroscopic molecu- 
lar vibrational data(see Ref.P| and references therein), 
the quartic anharmonic oscillator also serves as a ba- 
sic tool for checking various approximate and perturba- 
tive methods in quantum mechanics. Such an application 
appears in several recent field theoretical model studies 
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It is well known |15l [iG] that for the quartic an- 
harmonic oscillator the perturbation expansion diverges 
even for small couplings and becomes completely use- 
less for strong coupling. In view of this di verg ence 
of perturbation theory, we have adopted 0, 0| the 
general and very pow erful quasilinearization method 
(QLM) [11 13, m M El, which although iterative is 
not a perturbative method. In QLM the p-th order solu- 
tion of a nonlinear differential equation with N variables 
is obtained by first approximating the nonlinear aspect 



by a sequence of linear terms and then iteratively solving 
the associated linear equations. This iterative process 
converges to a solution without requiring the existence 
of a smallness parameter. Properties and applications of 
the quasilinearization method were reviewed recently in 

To apply the quasilinearization method, one first casts 
the Schrodinger equation into the nonlinear Riccati form 
and then solves that nonlinear equation by the QLM iter- 
ations. In a series of publications p^ ll8ll25ll26ll27| . we 
have shown that for a range of anharmonic and other 
physical potentials (with both weak and strong cou- 
plings), the QLM iterates display very fast quadratic 
convergence. Indeed, after just a few QLM iterations, 
energies and wave functions are obtained with extremely 
high accuracy, reaching 20 significant figures for the en- 
ergy of the sixth iterate even in the case of very large 
couphng constants. 

Although numerical solutions using either the QLM 
or direct numerical solution of the differential equations 
can be very accurate, it is important to also provide an- 
alytic solutions. Analytic solutions allow one to gauge 
the role of different potential parameters, and explore 
the influence of such variations on the properties of the 
quantum system under study. However, in contrast to 
the harmonic oscillator, the anharmonic oscillator can- 
not be solved analytically, and thus one usually has to 
resort to approximations. 

The goal of this paper is to obtain and test approxi- 
mate analytic solutions for the quartic and pure quartic 
oscillators using the explicit analytic equation for the first 
QLM iterate. We will show that both energies and wave 
functions will be represented by closed analytic expres- 
sions with the accuracy of the wave functions being be- 
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tween 0.1 and 1 percent for both small and large coupling 
constants. Various accurate analytic expressions for the 
energies have already appeared in the literature based on 
using convergent, strong coupling expansions generated 
by rearran gem ent of the usual divergent weak cou pling 
expansion |28l | or by some variational requirement |29l |. 
However, accurate analytic expressions representing wave 
functions have not hitherto been known. That result is 
provided here. 



II. MAIN FORMULAE 



The usual WKB substitution y{x) = ^j^^^ converts the 
Schrodinger equation to the nonlinear Riccati form 



dyjx) 
dx 



Here fc^(x) — 2[E — V{x)], where we use m = 1, fi. = 1 
units. The quasihnearization [U ill of Eq.Q 

leads to the recurrence differential equation 



dx 



+ {2yp^i{x))yp{x) = yl^i{x) - k'^{x), (3) 



where yp{x) is the subsequent p*'* QLM iterate, which 
have the same boundary condition as y{x) of Eq.(P)l. 
Note that Eq. (El is a linear equation of the form '^^^^^^ + 
f{x)yp{x) = q{x), with f{x) = 2yp^i{x) and q{x) = 
yl_i{x)-k^{x). 

Let us use Eqs.® to estimate the ground state wave 
function and energy of the quartic oscillator. Excited 
states will be considered elsewhere. 

The ground state wave function is nodeless and for an 
even potential should therefore be an even function. 
Its logarithmic derivative is necessarily odd, and there- 
fore the boundary condition obviously is y(0) = and 
correspondingly yp{0) = 0. 

A. Linear Initial Condition 

The zero iterate should be based on physical consider- 
ations. Let us consider first an initial guess yo{x) = —gx. 
This linear initial condition completely neglects the an- 
harmonic term containing A compared with the harmonic 
term and thus this initial guess is expected to be reason- 
able only for relatively small values of A. 

Solution of the first order linear differential Eq. Q with 
the above zero boundary condition at the origin can al- 
ways be found analytically. For p — 1 the solution is 



yi 



(x) = 2eS"'/ e-s'\g^s^ + Xs^~Ei)ds. (4) 
Jo 



Integration by parts, yields an expression for yi{x) that 
involves the error function Erf(x) = ^ e~* dt: 



yiix) 



+ e 



,5/2 



gx 



{-2x^ (2/ + 3X + 2gXx^) 
/^Erf(a;V^) [2g\g - 2E^) + i\] }. 
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The asymptotic expression Erf(a::) ~ (1 — for 

I a; I — > oo, indicates that yi{x) will be exponentially 
large for very large x unless the second term in Eq.Q 
is made zero. Correspondingly, invoking the condition 
2g^{g — 2E\) -|- 3A = yields the energy and the log- 
arithmic derivative in the first iteration: E\ — ^ and 
y\{x) = —ax — bx^, where a = g + b = ^. This 
leads to the first QLM iteration wave function ijj{x) — 



(2) C'exp 



gx_ 
2 



A 



This QLM result for the 



( 3^! 

energy coincides with the perturbative result, as well as 
with the result obtained by Friedberg, Lee and Zhao SO] 
who used their recently developed iterative method for 
solving the Schrodinger equation. 

The wave functions we obtained above obviously have 
incorrect asymptotic behavior. Also, the energies Ei cal- 
culated for different A, as displayed in Table 1, are far 
from being precise. Therefore, to improve the result one 
is tempted to go to the second QLM iteration, using yi {x) 
as an input. 

Eq.Q then yields the second iterate 



{[a^ + g^ + 2(ab + \)t^ (6) 



bH^]t^ - 2£:2}e""* "5* dt. 



Since y2{x)e^^°-^ ) approaches a constant when x 
goes to infinity, y2ix) and consequently the correspond- 
ing wave function grows exponentially at infinity, unless 
the integral in Eq.© equals zero when its upper limit 
equals infinity. This condition yields the following ex- 
pression for the energy E2 in the second iteration: 



E. 



2iab + X) t^ + bH^] t^e-"*^- 



5* dt 



-i^'dt 



(7) 

Values of i?2, for this initial linear form, are compared to 
exact values E calculated numerically in Table 1. It is 
seen that E2 approximates the exact E reasonably well 
only for small A, as we anticipated would result from 
using an initial linear condition. We now turn to another 
choice for the initial form. 



B. Quadratic Initial Condition 

To ensure a proper wave function asymptotically, one 
needs an adequate initial guess. Our second condition 
is based on the asymptotic behavior of the quartic po- 
tential. The zeroth iterate of the logarithmic derivative. 
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yo{x), is for example now obtained by taking an initial 
iterate of quadratic form yo{x) = — \/2A x"^. This choice 
for the first iterate automatically satisfies the asymptotic 
and X = boundary conditions, but is based on neglect- 
ing the harmonic term compare with the anharmonic one. 

With this quadratic initial condition, the solution 
yi (x) of Eq. (|3J) satisfying yi (0) = is then given by 



g§\/2A 



-4V2A 



[4A 



2 2 

9 s 



Note, that yi{x)e a 



- 2Ei] ds. 
(8) 

approaches a constant. Coo, 



as X goes to infinity, and consequently yi (x) grows expo- 
nentially at infinity unless the above limit Coo, is set to 
zero. The latter condition yields another expression for 
the energy based on the quadratic initial condition: 



El 



2Xt + 2Xr 



2V2X ^3 



dt 



Jo 



-at' 



dt 



(12) 



Another possible initial guess, which also accounts for 
the asymptotic behavior of both harmonic and anhar- 
monic parts of the potential is yo{x) — —y^g'^x^ + 2Xx^. 
This guess is easy to justify by assuming that y'{x) 
in Eq.© at large x becomes negligible compare with 
y'^{x) and therefore y'^{x) can be set equal k'^{x) = 
g^x^ + 2Xx^ — 2E K g^x^ -I- 2Xx'^ since in this expres- 
sion we can neglect E for sufficiently large x. 

The solution of Eq.Q using the initial condition 
yo{x) — —^g^x^ -I- 2Xx^, is only slightly more compli- 
cated than in when one uses the initial guess yo{x) = 
—gx — ^/2Xx^ , and is given by the expression 



El = 



roc 

Jo 



/2A 



2X s4 



ds 



Jo 



2A 



ds 



(9) 



The integrals in Eq. @ can be expressed in terms of the 
Euler Gamma function r(2;) |3lj . The final expression 
for the first iterate energy based on a quadratic initial 
condition reads 



El 



Ah 



-^ + x^33r - 
aVx V3 



(10) 



This result should be proper for large A when the sec- 

22 . 

ond term of the quartic potential ^-j — h Xx dominates 
over the harmonic term and thus a term containing g in 
initial guess yo{x) could be neglected. The above expres- 
sion for the energy is not expected to be suitable for small 
A. Indeed, one can see that for the A ^ 0, the energy in 
Ea. ((Tn|l diverges. 



C. Linear Plus Quadratic Initial Condition 

To obtain a result accommodating arbitrary values of 
both g and A, one has to start from an initial choice yo{x) 
that takes into account the asymptotic behavior of both 
the harmonic and the anharmonic parts of the potential. 
Let us choose, for example, yoix) = —gx— y/2Xx^ which 
is a linear combination of our two previous initial guesses. 
This yields 



£;i e^s* 3-* 



dt, 



g^ + g\/2Xt + 2Xt^ 



(11) 



with 



(s2+2Ax2) 2 

yi{x)=2e 3A 



(3^+2At^) 2 
3A 



[g^t^ + 2Xt^ - £^1] dt. 



(13) 

The condition to avoid exponential behavior of the log- 
arithmic derivative at infinity now yields the following 
expression for the ground state energy 



El 



roo,2 .2 



2Xt^ 



(32+2At2) 2 



dt 



Jo 



e 3A dt 



(14) 



For a pure quartic oscillator with g = 0, both Eas. l|ll|) 
and ifT^ reduce to 



2%/2Ax 

2/i(x) = 2e 3 
with El given by 

7.x I " p 

El = 



[2Xt^ - El] dt (15) 



nOO 2n/2A „fi 

Jg e 3 ' ds 



s^ds 1 1 r (1^ 

= A333. 



r(i) 



0.729011A3. 



(16) 



In view of Eq. (|16|l yi (x) of Eq. (|15|l can be expressed in 
terms of a special function,; namely, the Exponential In- 
e~''H-^'dt■. 



tegral EI^(^) = /i 



2Xx^ 



2 2V5a 3 



2Xx^'EIi 



2V2A 



Xhi-j4xEl2 

r(i) ^ 



2V2X 



}■ 



(17) 



This expression for the log derivative yields a first iter- 
ate QLM wave function, based on the initial condition 
yoix) = —\Jg'^x'^ -|- 2Ax^, that is the main result of our 
paper. 

The exact dependence of Ei on A for the pure quar- 
tic oscillator has the same form, but with a factor of 
0.667 986 259 before As so that the accuracy of 

the QLM prediction for the energy is about 9.1 per- 
cent. The WKB energy can be easily estimated and gives 
EwKB — 0.546267A3 , an accuracy of 18.2 percent. 
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TABLE I; Ground state energies E for the quartic oscillator 
potential U{x) = + Xx^ with g = 1. 



A 




EwKB 


^EwKB 


£(1) 




E^'^^wo = -gx 




£(3) 


AB*^' 


^(4) 


AB*"' 


yo = 


AB*^' 










yo = ~gx 


(%) 


2^^iteration 


(%) 


yo = -x^\/2A 


(%) 


yo = -gx - x^\/2A 


(%) 


-y^g'^x^ + 2\x'^ 


(%) 





1/2 


1/2 





1/2 





1/2 









1/2 





1/2 





0.1 


0.55915 


0.53328 


4.6 


0.575 


2.8 


0.55983 


0.1 


0.75658 


35 


0.56940 


1.8 


0.56149 


0.4 


0.3 


0.63799 


0.58466 


8.3 


0.725 


13.6 


0.64869 


1.7 


0.77799 


21.9 


0.64838 


1.6 


0.64705 


1.4 


0.5 


0.69618 


0.62538 


10.2 


0.875 


25.7 


0.72728 


4.4 


0.82319 


18.2 


0.70552 


1.3 


0.71126 


2.2 


1 


0.80377 


0.70420 


12.4 


1.25 


55.5 


0.91423 


13.7 


0.92313 


14.8 


0.81138 


0.95 


0.83090 


3.4 


2 


0.95157 


0.81667 


14.2 


2 


110 


1.2829 


34.8 


1.07257 


12.7 


0.95853 


0.73 


0.99577 


4.6 


10 


1.50497 


1.25412 


16.7 


8 




4.2628 


183 


1.6607 


10.3 


1.5259 


1.4 


1.61085 


7.0 


100 


3.13138 


2.57181 


17.9 


75.5 








3.4256 


9.4 


3.2564 


4.0 


3.40039 


8.5 


1000 


6.69422 


5.47955 


18.1 


750.5 








7.3095 


9.2 


7.1171 


6.3 


7.29744 


9.0 



III. RESULTS AND DISCUSSION 

The ground state energies for the quartic oscillator in 
the first QLM approximation for different initial guesses 
and for values of g = 1 and A between zero and one 
thousand and their comparison with the numerically cal- 
culated exact and WKB values are given in Table 1. One 
can see that the values computed using explicit equations 
(|12|l and H14|l for the QLM energy are significantly more 
accurate than the WKB values or than values obtained 
in the first and second QLM iterations with the initial 
guess yo{x) = —gx. They have a precision of 0.4 to 9 
percent for values of A varying between 0.1 and 1000, 
respectively. 

However, the main results of our work, are not the 
expressions for the energy. As mentioned in the intro- 
duction, such expressions were already given in different 
forms by others. Our major results are the analytic ex- 
pressions for the wave functions given by Eas. (|ll|l and 
Ea. (|13|l . which are based on using the first QLM iter- 
ate with the initial conditions yo = — \/2Ax^ and 
yo = — + 2Ax'', respectively. 

The graphs of the wave functions for the quartic os- 
cillator with g = 1 and for different A together with 
the correspondent exact and WKB wave functions, are 
presented in Figs. Q 13 and [Sj while Figs. EJ ^ and 
display the logarithm of the absolute value of the dif- 
ferences between the WKB or QLM wave functions and 
the exact solution for A being equal to 0.1, 1 and 10, 
respectively. The same graphs for the pure quartic os- 
cillator (g — 0) are presented in Figs. El El and in 
Figs. IHl El El respectively. One can see that in all the 
graphs the differences between the exact and QLM solu- 
tions are two to three orders of magnitude smaller than 
the differences between the exact and the WKB solutions 
and that the QLM wave functions expressed analytically 
by Eas. Hll|l . (|13|) have an accuracy of between 0.1 and 1 
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-x^sqrt(2>.) 


QLM; 


yo= 


-sqrt(g^x^+2Ax'') 


QLM; 
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FIG. 1: Comparison of the WKB, QLM and exact wave 
functions for the ground state of the quartic oscillator for 
ff = l,A = 0.1. 



percent. The order of magnitude better accuracy of the 
wave function compared to the poorer accuracy of the en- 
ergies is explained by the fact that the general theorems 
^2^2lLM, M, 2M for the QLM iterates show that 
the solutions converge quadratically with each iteration, 
while no such convergence theorem has been proven for 
the energy iterates. Note, that the dips in the Figures are 
artifacts of the logarithmic scale, since the logarithm of 
the absolute value of the difference of two solutions goes 
to — oo at points where the difference changes sign. The 
overall accuracy of the solution can therefore be inferred 
only at x values not too close to the dips. 
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FIG. 2: Logarithm of the differences of the WKB and QLM 
wave functions with exact wave function for the ground state 
of the quartic oscillator for g — 1,X — 0.1. 



FIG. 4: Logarithm of the differences of the WKB and QLM 
wave functions with exact wave function for the ground state 
of the quartic oscillator for g — 1, X — 1. 
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FIG. 3: Comparison of the WKB, QLM and exact wave 
functions for the ground state of the quartic oscillator for 
S = 1,A = 1. 



FIG. 5: Comparison of the WKB, QLM and exact wave 
functions for the ground state of the quartic oscillator for 
5 = 1,A = 10. 



IV. CONCLUSION 

We calculated analytically the ground state energy and 
wave function of the quartic and pure quartic oscillators 
by casting the Schrodinger equation into the nonlinear 
Riccati form , which is then solved in the first iteration of 
the quasilinearization method (QLM), which approaches 
the solution of the nonlinear differential equation by ap- 
proximating nonlinear terms with a sequence of linear 
ones and does not rely on the existence of a smallness 
parameter. Comparison of our results with exact nu- 
merical solutions and the WKB solutions shows that the 



explicit analytic expressions we obtain H12|l and (|14() for 
the ground state energy have a precision of only a few 
percent while the analytically expressed wave functions 
(|ll|l and (|13|l have an accuracy of between 0.1 and 1 
percent and are more accurate by two to three orders 
of magnitude than those given in the WKB approxima- 
tion. The QLM wave function in addition possess no 
unphysical turning point singularities which allows one 
to use these wave functions to make analytical estimates 
of the effects of variation of the oscillator parameters on 
the properties of systems described by quartic and pure 
quartic oscillators. 
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FIG. 6: Logarithm of the differences of the WKB and QLM 
wave functions with exact wave function for the ground state 
of the quartic oscillator for g — 1, X — 10. 



Extension of this approach to excited states and to 
other potentials is underway. 
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FIG. 7: Comparison of the WKB, QLM and exact wave func- 
tions for the ground state of the pure quartic oscillator for 
g = 0,\ = 0.1. 



FIG. 8: Logarithm of the differences of the WKB and QLM 
wave functions with exact wave function for the ground state 
of the pure quartic oscillator for g = Q,X = 0.1. 
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The next QLM iterations could be evaluated numeri- 
cally millilllllll. These further QLM iterates for 
the different anharmonic and other physical potentials 
with both strong and weak couplings also display very 
fast quadratic convergence so that the accuracy of ener- 
gies and wave functions obtained after a few iterations 
is extremely high, reaching 20 significant figures for the 
energy of the sixth iterate even in the case of very large 
coupHng constants. 



FIG. 9: Comparison of the WKB, QLM and exact wave func- 
tions for the ground state of the pure quartic oscillator for 
5 = 0,A = 1. 
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